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ABSTRACT 

We investigate the chemical properties of local galaxies within a cosmological frame- 
work in the hierarchical picture of galaxy formation. To this aim, we use a hierarchical 
semi-analytic model which includes the contribution from (i) low and intermediate 
mass stars, relevant producers of some important elements, such as C and N; (ii) type 
la Supernovae (SNe), for which a continuous delay-time distribution is assumed and 
which are important for the production of Fe and (iii) massive stars, dying as core- 
collapse type II SNe and which produce the a— elements. In this way, we can study 
abundances for a large set of chemical elements and in various galactic types, com- 
paring our predictions with available observations in the Milky Way (MW), in local 
dwarf galaxies and in local ellipticals. 

For Milky- Way-like galaxies, we can successfully reproduce the major observational 
constraints, i.e. the [O-Fe] vs [Fe/H] relation observed in disc stars and the stellar 
metallicity distribution (SMD). For dwarf galaxies, the stellar metallicity vs mass re- 
lation is reproduced by assuming that a substantial fraction of the heavy elements is 
lost through metal-enhanced outflows and a type la SN realization probability lower 
than the one of MW-like galaxies. The predicted abundance ratios for dwarf galax- 
ies are comparable to observations derived locally for dwarf spheroidals. The stellar 
metallicity distributions predicted for dwarf galaxies are in agreement with the lo- 
cal observations. We predict a substantial presence of extremely low metallicity stars 
[Fc/H]<-2.5, which have been recently observed in ultra-faint dwarf galaxies. In ellipti- 
cals, the observations indicate higher [a/Fe] values in larger galaxies. Several previous 
attempts to model the [a/Fe] vs a in ellipticals based on A— cold dark matter galaxy 
formation models predicted an anti-correlation between [a/Fe] vs a, indicating too 
much extended star formation histories in high-mass galaxies. Our results computed 
with a standard Salpeter initial mass function (IMF) indicate a flat [a/Fe] vs relation. 
However, we suggest a possible solution to this problem and show how, by assuming 
a star formation-dependent IMF with a slope x = 1.35 in systems with star formation 
rates < 100Af Q /?/r and slightly flatter (i.e. with x = 1) in object with stronger star for- 
mation, the observed correlation between [a/Fe] and a can be accounted for on a large 
velocity dispersion range. Fundamental roles are played also by interaction-triggered 
starbursts and AGN feedback. Finally, a SF-dependent IMF seems necessary also to 
reproduce the stellar metallicity-er relation observed in local early-type galaxies. 

Key words: Galaxies: formation and evolution; Galaxies: abundances. 



1 INTRODUCTION 

In the last few years, the power of high-resolution spectro- 
graphs has made available a large amount of stellar and 
interstellar abundances for various chemical species. Nowa- 
days, large catalogues of stellar abundances are accessible for 
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the Milky Way galaxy and for the dwarf galaxies of the Lo- 
cal Group. Gas-phase abundances have been derived up to 
high-redshift (2 > 3), both by analysing emission lines and 
absorption lines, present in the spectra of distant Quasars. 
This large set of data provides us with valuable constraints 
for galactic chemical evolution studies. Of particular interest 
is the study of chemical evolution by means of "ab-initio" 
galaxy formation models. 
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The study of galactic chemical evolution within a cosmolog- 
ical framework is of particular importance, since it provides 
us with crucial information on the star formation history 
(SFH) of galaxies, on the ages of the stellar populations and 
on the gas accretion and outflow histories. Chemical evolu- 
tion offers a way to constrain the main parameters driving all 
of these processes, which are of primary interest for galaxy 
formation theories. 

So far, a few theoretical studies of hierarchical galaxy forma- 
tion have payed attention to the problem of galactic chemical 
evolution. In a pioneeristic paper, Thomas (1999) considered 
star formation histories from a hierarchical semi-analytic 
galaxy formation model, taking into account for the first 
time the rates of type la SNe. This attempt allowed him to 
compute the evolution of the abundance ratio between ele- 
ments formed on different timescales, i.e. O, synthesised by 
type II Supernovae (SNe) on timescales less than 0.03 Gyr, 
and Fc, produced mainly by type la SNe on timescales rang- 
ing from 0.03 Gyr up to one Hubble time. Other important 
works on chemical evolution in the hierarchical framework 
are the ones by Nagashima et al. (2005), and Nagashima 
& Okamoto (2006), where the chemical evolution of vari- 
ous elements produced by type la and type II SNe is com- 
puted. However, none of the works quoted above takes into 
account the role of low and intermediate mass stars (LIMS 
hereinafter), i.e. all the stars with masses 0.8 < m/M Q < 8. 
These stars are of fundamental importance in chemical evo- 
lution studies since they incorporate an important fraction 
of the metals and during the last stages of their evolution, i.e. 
during the planetary nebula phase, they may restore vari- 
able amounts of non-processed heavy elements incorporated 
at their birth and produced by stars of previous generations, 
giving a significant contribution to the heavy element pol- 
lution of the interstellar medium (ISM). Furthermore, low 
and intermediate mass stars are relevant producers of C and 
N (Van den Hoeck and Groenewegen 1997), whose abun- 
dances can not be properly assessed by any chemical evolu- 
tion model which does not take into account in detail the 
evolution of LIMS. In this paper, with a hierarchical semi- 
analytical model (SAM) we study the chemical evolution of 
these heavy elements, taking into account also the contribu- 
tion from LIMS. Our main aim is providing chemical evolu- 
tion predictions computed by means of star formation (SF) 
and gas accretion histories derived from an ab-initio galaxy 
formation model. For the first time, by means of a SAM 
for galaxy formation, we study the evolution of the abun- 
dances for a large set of chemical elements produced by stars 
of various masses, ending their lives on various timescales. 
This paper is mainly focused on the chemical abundances in 
local galaxies. We compare our predictions to abundances 
observed in various galactic environments, such as in the 
Milky- Way disc, local dwarf galaxies and local early-type 
galaxies. Future papers will be devoted to the study of chem- 
ical abundances at high redshift. This paper is organised as 
follows. In Sect. 2, we present our theoretical instruments, 
i.e. the hierarchical semi-analytical galaxy formation model 
and the main chemical evolution equations. In Sect. 3, we 
present our results. Finally, in Sect. 4 we discuss the main 
implications of our results and draw some conclusions. 



2 THE MODEL 

In this paper, we start from the star formation histories of 
galaxies computed by means of a semi-analytic model of 
galaxy formation. For each star formation history, we com- 
pute the chemical evolution a posteriori by means of detailed 
chemical evolution equations. In this section, we briefly de- 
scribe the SAM used in this work and the methods used to 
compute the evolution of the chemical abundances. 



2.1 The galaxy formation model 

We derive the star formation histories of galaxies and of 
their progenitors from the semi-analytic model by Menci et 
al. (2006, 2008); we recall here its key features. Galaxy for- 
mation and evolution is driven by the collapse and growth 
of dark matter (DM) halos, which originate by gravitational 
instability of overdense regions in the primordial DM den- 
sity field. This is taken to be a random, Gaussian density 
field with Cold Dark Matter (CDM) power spectrum within 
the "concordance cosmology" (Spergel et al. 2006) for which 
we adopt round parameters Qa = 0.7, Q,o = 0.3, baryonic 
density fij, = 0.04 and Hubble constant (in units of 100 
km/s/Mpc) h = 0.7. The normalisation of the spectrum is 
taken to be as = 0.9 in terms of the variance of the field 
smoothed over regions of 8 h' 1 Mpc. 

As cosmic time increases, larger and larger regions of the 
density field collapse, and eventually lead to the formation of 
groups and clusters of galaxies, which grow by merging with 
mass and redshift dependent rates provided by the Extended 
Press & Schechter formalism (see Bond et al. 1991; Lacey 
& Cole 1993). The clumps included into larger DM halos 
may survive as satellites, or merge to form larger galaxies 
due to binary aggregations, or coalesce into the central dom- 
inant galaxy due to dynamical friction; these processes take 
place on timescales that grow longer over cosmic time, so the 
number of satellite galaxies increases as the DM host halos 
grow from groups to clusters. All the above processes are 
implemented in our model as described in detail in Menci et 
al. (2005, 2006), based on canonical prescriptions of semi- 
analytic modeling. 



2.1.1 Star Formation and SN feedback 

The radiative gas cooling, the ensuing star formation and 
the Supernova events with the associated feedback occur- 
ring in the growing DM halos (with mass M and circular 
velocity v) are described in Menci et al. (2005). The cooled 
gas with mass M g settles into a rotationally supported disk 
with radius r d (typically ranging from 1 to 5 kpc), rotation 
velocity v d and dynamical time td = r d /v d . The gas grad- 
ually condenses into stars and the stellar feedback returns 
part of the cooled gas to the hot gas phase with mass 
at the virial temperature of the halo. As for the star forma- 
tion, we assume the canonical Schmidt form tp — M g /(qT d ), 
where r d = r d /v d and q is fixed by the Kennicutt (1998) 
law. At each time step, the mass Ami, returned from the 
cold gas content of the disk to the hot gas phase due to Su- 
pernovae (SNe) activity is estimated from canonical energy 
balance arguments (Kauffman 1996, Kauffmann & Chariot 
1998; see also Dekel & Birnboim 2006) as 
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Am h = EsncoVo Am*/« 2 



(1) 



where Am* is the mass of stars formed in the timestep, 
ij « 3 — 5 ■ 10~ 3 /M Q is the number of SNe per unit solar 
mass (depending on the assumed IMF), Esn = 10 51 ergs is 
the energy of ejecta of each SN, and v c is the circular veloc- 
ity of the galactic halo; e = 0.01 — 0.5 is the efficiency of the 
energy transfer to the cold interstellar gas. The above mass 
Am/, is made available for cooling at the next timestep. The 
model free parameters q = 30 and eo = 0.1 are chosen as to 
match the local B-band luminosity function and the Tully- 
Fisher relation adopting a Salpeter IMF. 
Note that our simple modelling for SNae feedback provides 
a very good fit to the observed correlations of outflow ve- 
locities with galactic properties (like the circular velocity or 
the star formation rate). Although an accurate estimate of 
the outflow velocity would require a detailed treatment of 
the gas kinematics (including also the physics of OB asso- 
ciations and the dynamics of of superbubbles, see Ferrara 
et al. 2000; Veilleux et al. 2005) and is beyond the aim of 
the present paper, its asymptotic value V ou fiow can be esti- 
mated from the energetic balance between the total energy 
deposition rate by SNe and the rate of kynethic energy loss 
as (Veilleux et al. 2005) Ami V^ OI = Esn eo Vo tp- A com- 
parison with our assumed expression for the expelled gas 
mass Arrih given above implies V m ax ~ v c (or equivalently 
a mass loading factor rhh/ip cx v~ 2 ), a relation which is 
in very good agreement with the relation observed in local 
starbursts (Martin 2005; Rupke et al. 2005; Veilleux et al. 
2005). 



2.1.2 AGN feedback 

The model also includes a treatment of the growth of super- 
massive black holes at the centre of galaxies by interaction- 
triggered inflow of cold gas, following the physical model of 
Cavaliere & Vittorini (2000). Our SAM includes a detailed 
treatment of feedback from Active Galactic Nuclei (AGN), 
which acts only during the active AGN phases of each galaxy 
(hence for a minor fraction ~ 10~ 2 of the galaxy lifetime). 
This is assumed to stem from the fast winds with veloc- 
ity up to 10 _1 c observed in the central regions of AGNs 
(Weymann 1981; Turnshek et al. 1988; Risaliti et al. 2005); 
these are usually though to originate from the acceleration 
of disk outflows due to the AGN radiation field (Proga 2007 
and references therein, Begelman 2003). These supersonic 
outflows compress the gas into a blast wave terminated by 
a leading shock front, which moves outwards with a lower 
but still supersonic speed and sweeps out the surrounding 
medium. Eventually, this is expelled from the galaxy. Quan- 
titatively, the energy injected into the galactic gas in such 
inner regions is taken to be proportional to the energy ra- 
diated by the AGN, AE — eagn rjc 2 Am acc . the value of 
the energy feedback efficiency for coupling with the sur- 
rounding gas is taken as cagn = 5 • 10~ 2 , consistent with 
the values required to match the X-ray properties of the 
ICM in clusters of galaxies (see Cavaliere, Lapi & Menci 
2002). This is also consistent with the observations of wind 
speeds up to v w ~ 0.1c in the central regions, that yield 
cagn ~ v w /2c ?s 0.05 by momentum conservation between 
photons and particles (see Chartas et al. 2002, Pounds et 
al. 2003); this value has been also adopted in a number of 



simulations (e.g., Di Matteo, Springel & Hernquist 2005) 
and semi-analytic models of galaxy formation (e.g., Menci 
et al. 2006). The transport of the above energy in the galaxy 
and its effect on the distribution of the galactic gas is com- 
puted in detail by describing the expansion of the blast wave 
solving the corresponding hydrodynamical equations; these 
include the effects not only of initial density gradient, but 
also those of upstream pressure and DM gravity (Lapi, Cav- 
aliere, Menci 2005). The solutions show in detail how the 
perturbed gas is confined to an expanding shell bounded by 
an outer shock at the radius R s {t) which sweeps out the gas 
surrounding the AGN. From the resulting shock expansion 
law Rs(t) the amount of expelled gas can be computed (for 
different AGN luminosities and galactic properties), as well 
as the final galactic gas distribution, as described in detail 
in Menci et al. (2008) We refer to Menci et al. (2005-2008) 
for details and for the comparison of the model results with 
observations concerning the cosmological evolution of both 
the galaxy and the AGN population. 



2.1.3 Starbursts triggered by galaxy interactions 

An additional channel for star formation implemented in 
the model is provided by interaction-driven starbursts, trig- 
gered not only by merging but also by fly-by events be- 
tween galaxies; such a star formation mode provides an 
important contribution to the early formation of stars in 
massive galaxies, as described in detail in Menci et al. 
(2004, 2005). The galaxy interaction rate is given by t^ 1 = 
nr S(r t , v, V re i) Vrei, where tit is the number of galaxies 
hosted in a given dark matter halo, E is the cross section 
for grazing (i.e. at distances closer than the galaxy tidal ra- 
dius) encounters, computed as in Menci et al. (2003) from 
the orbital parameters (the impact parameter b, the ra- 
dius of the host dark matter halo) for each galaxy in our 
Monte Carlo simulations; the relative velocity V re i is com- 
puted from the velocity dispersion of the dark matter halo 
hosting the interacting galaxies. With a probability given 
by the above equation in each time step, a galaxy is con- 
sidered to be in a starburst phase. In such a case, the 
cold gas fraction converted into stars during the burst is 



face- f- = ICT 1 ^ r -f^y The fraction of cold gas 

accreted by the black hole in an interaction event is here 
computed in terms the variation Aj of the specific angular 
momentum j ~ Gm/vd, induced by the grazing encounter. 
Here m' is the mass of the partner galaxy in the interaction, 
and the average runs over the probability of finding such a 
galaxy in the same halo where the galaxy with mass m is 
located. The duration of the burst is given by the interac- 
tion timescale Vd/V re i. With such a description, a fraction of 
gas « 10% is impulsively converted into stars during minor 
bursts with m'/m ~ 0.1, while in the rarer major bursts 
with m'/m w 1 the fraction of gas converted into stars can 
approach 100%. 



2.2 Galactic chemical evolution 

In this paper, we model chemical evolution following the 
approach of Matteucci & Greggio (1986). We follow the time 
evolution of the abundances for the following set of chemical 
elements: H, He, C, N, O, Mg, Si, Fe. This is one of the major 
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novelties of this work, since so far such a large set of chemical 
elements has never been considered in semi analytic models 
for galaxy formation. 

For any chemical element i, the variation of its mass in the 
interstellar gas per unit time is: 

" !>U = -Xi(t)^{t) + R.LIMS,i{t) + Rla,i(t) + 



dt 



Rii,i(t)+Ii(t) + Oi(t) 



(2) 



where Xi is the mass fraction in the gas for the element i 
at the time t, ip(t) is the star formation rate (SFR) of the 
selected galaxy. The terms RLiMS,i(t), Ri a ,i(t) and Rn,i(t) 
are the rates of production of the element i from stars of 
low and intermediate mass, type la SNe and type II SNe, 
respectively. 

Finally, the terms Ii(t) and Oi(t) take into account possible 
mass increments and outflows of mass in the form of the 
element i, respectively. The mass increment may occur by 
means of three processes: (i) cooling of hot gas; (ii) by means 
of cool gas accretion or (iii) merging with other galaxies, 
whereas an outflow can occur owing to a galactic wind. In 
the remainder of the paper, when we will use the expression 
"Infall" , we will intend actually the increment of gas, which 
can occur by means of any of the three processes mentioned 
above. The methods used to compute the accretion and out- 
flow histories of the galaxies will be described in Sect. 2.3. 
Stars of low and intermediate mass are those with initial 
masses 0.8 < tti/Mq < 8. Their production rate for the 
element i is given by: 

fM(t) 

RhiMS,i{t)= I ip(t - T m )iriLiMS,i(m)<p(m)dm (3) 



-s 

JO.! 



where M(t) is the turnoff mass at the time t, r m is the life- 
time of the star of mass m, <j)(m) is the initial mass function 
(IMF) and m,LiMS,i{m) is the total mass in the form of the 
the element i that a low or intermediate mass star restores 
into the ISM. 

Unless otherwise stated, for the initial mass function, we as- 
sume a Salpeter (1995) law, given by <j)(m) = C ■ m~ ( - 1+x ) 
with x — 1.35. The constant C is determined by the condi- 
tion 

I m(f>(m)dm — l (4) 
Jo.i 

i.e. from the normalisation to 1 of the IMF by mass. The 
quantities iriLiMS,i{m) are taken from Van den Hoeck& 
Groenewegen (1997), who computed theoretical yields for 
various elements and for stars of low and intermediate mass 
as a function of the initial metallicity Z. The yield Pi(m) is 
defined as the fraction of the initial mass that a star trans- 
forms into the chemical element i and restores into the ISM. 
The quantity rriLiMS,i(m) is computed from the yield ac- 
cording to 



m LIM .s,i(m) = Xi(t - T m )m ej (m) + mpi(m) 



(5) 



where Xi(t — r m ) and m ej (m) are the original abundance 
of the element i computed at the time t — r m , i.e. when the 
star of mass m formed, and the total ejected mass for a star 
of initial mass m, respectively. 

For the element i, the production rate by type la SNe is 
Ria,i{t) — mi at iRia{t). Ri a {t) is the type la supernova (SN) 
rate at the time t, computed according to: 



Ria(t) = k a Ai 



rnin(t,T x ) 



i/j(t - r)DTD(r)dr, 



(6) 



(Greggio 2005; Matteucci et al. 2006). The quantity k a is 
the number of stars per unit mass in a stellar generation, 
given by: 



k a 



4>(m)dm 



(7) 



where m; = O.IM© and Ml = WOMq, whereas Ai a is the 
realization probability for type la SNe and represents the 
fraction of binary systems which may give rise to type la 
SNe, assumed to be constant in time. This quantity is uncer- 
tain (Maoz 2008) and is likely to depend on the environment 
(Matteucci et al. 2006). In chemical evolution studies, Ai a 
is basically treated as a free parameter, tuned in order to re- 
produce the type la SN rate in local galaxies (see Matteucci 
et al. 2006; Calura & Matteucci 2006). We will test the im- 
pact of this parameter on some of the results described in 
this paper. 

The function DTD(t) is the delay time distribution. We as- 
sume that Type la SNe originate from the explosion of a 
C/O white dwarf (WD) in a close binary system, where the 
companion is either a Red Giant or a Main-Sequence star. 
This is the single-degenerate (SD) model of Whelan & Iben 
(1973). The DTD(t) is from Matteucci et al. (2006 see also 
Matteucci & Recchi 2001, Valiante et al. 2009). 
mia t i is the mass that a type la SN synthesises and ejects 
into the ISM in the form of the element i. These quantities 
are taken from Iwamoto et al. (1999). 

We assume that all stars with mass m > 8Mq explode as 
type II SNe. The rate of production of the element i from 
type II SNe is: 

(•100 

Rn,i(t) = / ip(t - Tm)mii l i(m)cj>(m)dm 



f 

J8 



mi i ti (m)(j)(m)dm 



(8) 



the above simplification means to assume the instantaneous 
recycling approximation to treat chemical enrichment from 
massive stars, i.e. to assume that all stars with masses 
m > 8Mq have very small lifetimes. This assumption is 
motivated by the fact that the timestep we use to compute 
chemical evolution is of ~ 0.015 Gyr, i.e. half the lifetime 
of a massive star of 8Mq. With such a timestep, the effects 
of the lifetimes of single massive stars of various masses can 
not be appreciated. 

Chemical enrichment from type II SNe is computed by 
means of the quantity m//,i 



mu t i(m) — mii,new,i + rriej • Xi{t) 
where 



mii., 



■■ pi(m)m 



(9) 



(10) 



The quantity pi (m) is the fraction of the mass m tranformed 
into the element i and is commonly known as the yield. The 
second term on the right side of eq. 9 is the mass in the form 
of the element i already present when the star formed: m ej 
is the total mass ejected by a star of mass m and Xi is the 
interstellar mass fraction of the element i at the time t. For 
massive stars, we adopt the metallicity- dependent yields of 
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Woosley & Weaver (1995). 



2.3 Determining the past star formation, 
accretion and the outflow history 

For any galaxy at redshift z ~ 0, it is possible to recon- 
struct its star formation rate ip(t) and its accretion and out- 
flow histories on the basis of its past merging history. We 
have already described the possible gas accretion processes 
in Sect. 2.2. The outflows may originate owing to the con- 
stribution of SN explosions and of the AGN feedback. As 
described in Section 2.1.1 and 2.1.2, both mechanisms in- 
ject energy into the ISM and cause gas to move from the 
cold phase to the hot phase. The hot gas may eventually be 
ejected if its thermal energy is larger than its binding en- 
ergy, due to the gravitational well of both the dark matter 
and the baryons. For a particular galaxy, its star formation 
rate at a past time t is given by 

^(t) = E^i(*) (11) 

where ipj (t) is the SFR of the j-th progenitor of the galaxy 
at the time t. This formula allows us to take into account the 
complex merging history of galaxies, playing an important 
role in chemical evolution. In a similar way, at any past 
time t, for any galaxy we are able to determine the total 
cold gas mass M g (t). If A gas (t) = M g (t + dt) - M g (t) - 
~EiR*,i ■ dt + ip{t)dt is the difference in the total gas mass 
at two following timesteps minus the total mass restored by 
the stars in dt plus the mass ending in stars, we use this 
quantity to compute the total "increment" and outflow rate 
at any timestep. If A ga s{t) / Mq > then we assume that at 
the time t 

I{t) = A gas /dt ; O(t) = (12) 
otherwise, 

l(t) = ; 0(t) = A gas (t)/dt (13) 

By means of these quantities, we are able to solve the chemi- 
cal evolution equations described in Sect. 2.2. For any chem- 
ical element i, Ii(t) = Xi(t)I(t) in case of merging or 
Ii(t) = X itP (t)I(t) in case of infall of pristine gas, where 
Xi iV is the primordial mass fraction for the element i. 

In The SAM, the metal-enriched gas expelled from 
satellite galaxies residing in a common DM halo is diluted by 
the hot gas of the halo, whose amount is much larger, gener- 
ally by a factor 100, than the amount of the interstellar cold 
gas. For this reason, the metallicity of the gas which may 
eventually fall back from the hot halo onto a galaxy, which in 
general is the central galaxy of the halo, may be considered 
negligible. This allows us to assume that the accreted gas has 
a primordial composition, i.e. Xh, p = 0.75, XHe, P = 0.25 
and for any other element, we assume Xi iP — 0. On the 
other hand, for any element, unless otherwise stated, we as- 
sume Oi(t) — Xi(t)0(t), i.e. the chemical composition of 
the outflow is the same as the one of the cold ISM. 



3 RESULTS 

In this section, we describe the results obtained for various 
types of galactic systems. These results concern mainly lo- 



cal galaxies and are compared to observational constraints 
obtained in the Milky Way, in local dwarf galaxies and in 
local ellipticals. A study of the chemical evolution of distant 
galaxies is the subject for future work. 

At z ~ 0, from the catalogue of galaxies obtained by means 
of the SAM of Menci et al. (2008), we select galaxies on 
the basis of various criteria. The selection criteria are simi- 
lar to those adopted for the observations, and are based on 
the combined use of various available quantities, such as the 
colour, star formation rate or the stellar mass. In the follow- 
ing sections, we describe how we select galaxies of various 
types and our comparison with observational data from the 
literature. 



3.1 Chemical abundances in Milky- Way like 
galaxies 

The first important step of our study is to test our pre- 
dictions for Milky- Way (hereinafter MW) like galaxies. Re- 
cently, De Rossi et al. (2009) published a detailed study 
of MW-like galaxies extracted from the Millennium Simu- 
lation catalogue. Following De Rossi et al. (2009), we de- 
fine MW-like galaxies as those with circular velocity V c in 
the range 200 < V c /km s^ 1 < 240 and cold gas fractions 
f c = M g /(M 3 + M.) in the range 0.1 < f c < 0.3. 
By means of these criteria, for the local number density of 
MW-like galaxies, we obtain ~ 1 • 10~ 4 Mpc~ 3 , consistent 
with De Rossi et al. (2009). 

In Fig. 1, we show the star formation histories of a set of 
9 MW-like galaxies drawn from our sample. In Fig. 2, wc 
show the cumulative stellar and gas masses as a function of 
time for these galaxies, whereas in Fig. 3 we show the time 
evolution of the gas accretion rates. 

The predicted abundance ratios between any chemical el- 
ement and Fe are sensitive to the parameters we have as- 
sumed to describe the type la SN rate. In particular, having 
assumed in this case a Salpeter IMF constant in time, the 
most relevant parameter in this study is the type la real- 
ization probability Aj a . In the literature, this quantity is 
not well-constrained, but it is generally considered as a free 
parameter tuned in order to reproduce the local type la 
SN rate (Matteucci et al. 2006; Calura & Mattcucci 2006). 
In Fig. 4, we show the present-time type la SN rate, ex- 
pressed in SNuMt, predicted for our selected MW-like galax- 
ies and compared to an observational range for the Milky- 
Way global la SN rate. The observational range is computed 
by assuming a SN rate of 0.3-0.4 1/cen (van den Bergh 
& Tamman 1991) and a present-day disc stellar mass of 
~5x 10 10 M Q (Mera et al. 1998). This figure shows that 
the assumption of a la SN frequency Ai a = 0.002 provides 
present-day type la SN rates in reasonable agreement with 
the observations (upper panel of Fig. 4) whereas by assum- 
ing Ai a = 0.004 we severely overestimate the observed type 
la SN rate (lower panel of Fig. 4). From this moment on, 
Ai a = 0.002 be our reference value for the type la SN real- 
ization probability. 

The chemical abundances obtained for the galaxies selected 
according to the criteria described at the beginning of this 



t 1 SNuM=l SN century- 1 10- 10 M® 
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Figure 1. Star Formation histories of a few MW-like galaxies selected according the criteria described in Sect. 3.1. 
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Figure 2. Evolution of the cumulative stellar masses (solid lines) and gas masses (dotted lines) for a few MW-like galaxies selected 
according the criteria described in Sect. 3.1. 
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Figure 3. Time evolution of the increment rate I(t) for a few MW-like galaxies selected according the criteria described in Sect. 3.1. 
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Figure 4. Predicted type la SN rate, expressed in SNuM (1 SNuM = 1 SN century -1 10~ 10 Mq) for selected MW-like galaxies compared 
to an estimate of the type la SN rate for the Milky Way (black solid lines). The colour code, shown by the bar at the top of the Figure, 
represents the predicted number of galaxies with a given SNuM and a stellar mass M 4 , normalised to the total number of galaxies with 
that stellar mass. In the upper and lower panels, the SN rate is computed by assuming Af a = 0.002 and Aj a = 0.004, respectively. 
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Figure 5. Predicted distribution of the abundance ratios [X/Fe] vs [Fe/H] for selected MW-likc galaxies for several chemical elements 
compared with observations from various authors. The colour code, shown by the bar at the top of the Figure, represents the predicted 
number of stellar populations belonging to MW-like galaxies with a given abundance ratio at metallicity [Fe /H] , normalised to the total 
number of stellar populations at that metallicity. Observational data: asterisks: Cayrel et al. (2004); plus signs: Carbon et al. (1987); 
open squares: Spite et al. (2005), open triangles: Laird (1985a; 1985b), the crosses: Tomkin et al. (1995), open diamonds: Clegg et al. 
(1981), solid circles: Israelian et al. (2004), solid diamonds: Frangois et al. (2004), finally the solid squares are from a compilation of data 
by Cescutti (2008). 



section are now compared to the abundances observed in 
Milky Way stars. 

In Fig. 5, we show the predicted abundance ratios vs 
metallicity for our sample of MW-like galaxies, compared 
to a set of observations of abundance ratios in MW stars 
for various elements. The predictions are represented by 
the colour code at the top of Fig. 5, expressing the pre- 
dicted number of stellar populations belonging to galaxies 
with given abundance ratio at metallicity [Fe/H], normalised 
to the total number of stellar populations at that metallic- 
ity. The symbols are the observational values from several 
sources (for further details, see caption of Fig. 5). In Fig. 5, 
the observational data are single-star abundances, i.e. each 
point represents the abundance ratios observed in a single 
local star. On the other hand, in our model we cannot resolve 
single stars, hence we consider the average abundances in a 
stellar population and we compare them to observations. For 
each selected MW-like galaxy, in the time interval dt a stellar 
mass dM* will be formed, characterized by a given metallic- 
ity (i.e. [Fe/H]) and particular abundance ratios. In Fig. 5, 
the coloured regions represent the predicted [X/Fe] - [Fe/H] 
distributions computed considering each stellar population 
of mass dM* born in each selected MW-like galaxy. This is 
very similar to what is done in chemical evolution models 
for the solar neighbourhood (see Matteucci 2001), barring 
the fact that in that case single tracks, representing the pre- 
dictions for one single star formation history, are plotted in 



the [X/Fe]-[Fe/H] diagrams. In our case, instead of having 
one single evolutionary track, we have various model tracks, 
each one drawn from a particular star formation history for 
a MW-like galaxy. All these tracks populate the coloured 
areas of Fig. 5 . 

It is important to note that we cannot explore the metal- 
licity region [Fe/H]< —2.5. To study chemical evolution at 
metallicities lower than this value, it would be necessary to 
consider finite stellar lifetimes also for massive stars. How- 
ever, the metallicity range —2.5 < [Fe/H] < encompasses 
the Fe abundances of all thin and thick disc stars, hence 
it is useful to investigate the chemical evolution of MW- 
like galaxies during their major disc phases. In Fig. 5, in the 
[O/Fe] vs [Fe/H] plot the predictions indicate the presence of 
two main populations of stars. The first population has the 
[O/Fe] ratio in anti-correlation with the [Fe/H] and overlaps 
with the observation of local stars. Beside this, Fig. 5 shows 
the presence of a horizontal population of stars with metal- 
licity -2. < [Fe/H] < -0.2 and constant [O/Fe]- 0.3. 

A similar behaviour is visible also in the [C/Fe] — [Fe/H], 
[Mg/Fe] — [Fe/H], and [Si/Fe]-[Fe/H] plots. This is due to 
substantial late increment episodes, very frequent in selected 
MW-like galaxies (see Fig. 3). The late accretion of pristine 
gas has the effect of decreasing the metallicity [Fe/H] of each 
stellar population born immediately after the infall event, 
having little effect on the [O/Fe] ratio, since both O and 
Fe are diluted by the same amount. Note that none of the 
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Figure 6. From top-left, clockwise: [a/Fe] versus time, [Fe/H] versus time, Infall rate versus time and [a/Fe] versus [Fe/H] for one 
the selected MW-like galaxies discussed in Section 3.1 
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Figure 7. Star Formation histories of MW-like galaxies whose mass assembled at z = 2 is at least 75% of their present mass. 
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Figure 8. Time evolution of the increment rate I(t) for MW-like galaxies whose assembled mass at z = 2 is at least 75% of their present 
mass. 
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Figure 9. Predicted distribution of the abundance ratios [X/Fe] vs [Fe/H] for MW-like galaxies whose mass assembled at z = 2 is at least 
75% of their present mass, for several chemical elements compared with observations from various authors. The colour code, shown by 
the bar at the top of the Figure, represents the predicted number of stellar populations belonging to MW-like galaxies born with a given 
abundance ratio and with metallicity [Fe/H], normalised to the total number of stellar populations with that metallicity. Observational 
data as in Fig. 5 
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Figure 10. Effects of varying the stellar yields and the IMF on the predicted [C/Fc] vs [Fe/H] for selected MW-like galaxies. We show 
our results obtained assuming a Salpeter IMF and the standard stellar yields described in Sect. 2. 2 (upper panel), a Kroupa et al. (1993) 
IMF and standard yields (middle panel), a Kroupa et al. (1993) IMF and the C yields of low and intermediate stars increased by a factor 
2. Colour code and observational data as in Fig. 5. 
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Figure 11. Predicted [O/Fe] vs [Fe/H] for Milky Way galaxies selected on the basis of their gas accretion history (upper panel) and for 
the ones presenting a ratio between the stellar mass in the disc and total stellar mass M, d/M*,tot = 0.7 — 0.8 (lower panel). Colour code 
and observational data as in Fig. 5. 
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Figure 12. Predicted cumulative stellar metallicity distribution of selected Milky-Way like galaxies (solid line) compared to local 
observations from various authors. The model predictions and the observational data are normalised by imposing that the subtended 
areas have all the same value, equal to unity. Observations are from Wyse & Gilmore (1995) (crosses); Rocha-Pinto & Maciel (1996) 
(triangles) and Jorgcnsen (2000) (diamonds) 



observational data in Fig. 5 shows this peculiar behaviour. 
To better understand this point, in Fig. 6 we show for one 
of the selected MW-like galaxies the time evolution of the 
average interstellar [a/Fe] and [Fe/H], of the total infall rate 
and the [a/Fe] vs [Fe/H] diagram. From the [Fe/H] vs time 
and infall vs time plots, one can clearly see that in corre- 
spondence of two major gas accretion episodes at 7.5 Gyr 
and ^11 Gyr, the [Fe/H] drops abruptly whereas the [a/Fe] 
remains unchanged. In the [a/Fe] - [Fe/H] diagram, the late 
accretion episodes and the consequent decrease of interstel- 
lar [Fe/H] are visible in the 'turn-off ' of the track at [Fe/H] 
~ 0, where the curve bends back towards lower [Fe/H], at 
constant [a/Fe] ~ -0.5. The decrease of [Fe/H] at constant 
[a/Fe] stops at [Fe/H] ~ -1.7, where the track bends for- 
ward, moving towards higher [Fe/H] values. 
To investigate this point, we perform a further selection 
on the sample of MW-like galaxies. The basic reason of 
the presence of the horizontal population of stars in the 
[0/Fe]-[Fe/H] is late accretion episodes. Now we focus on 
a given redshift, in order to understand which fraction of 
the present-day mass must be assembled in galaxies pre- 
senting chemical features the most similar to those of the 
MW as possible. This test may be useful to give constraints 
on the gas accretion history of the MW galaxy. We may 
focus on z — 2, which may be used a lower limit to the 
highest redshift of major merging for MW-like galaxies (see 
e.g. Colavitti et al. 2008). It is necessary that at least 75 
% of the present-day mass must be assembled at z = 2, in 
order to have the horizontal population vanished. In Fig. 7, 
we show the SF histories of these galaxies selected on the 



basis of their accretion history. The accretion histories of 
these galaxies are shown in Fig. 8. Basically, these galaxies 
do not experience major accretion episodes at epochs after 
~ 6 Gyr, at variance with most of the systems represented in 
Fig. 3, whereas their star formation histories show low star 
formation rates at late times. In Fig. 9, we show the the pre- 
dicted abundance ratios vs metallicity for all these galaxies. 
From Fig. 9, we see that in most of the cases, the observed 
behaviour of the abundance ratios measured in Milky Way 
stars is satisfactorily reproduced. 

Our results suggest that the abundance ratios observed in 
the Milky Way stars exclude that our Galaxy experienced 
major pristine gas accretion episodes within the last few 
Gyrs, either due to merging with gas-rich, low-metallicity 
dwarf galaxies, or due to infall of pristine gas. 
Our results are in agreement with a large set of previous pa- 
pers, indicating that the merging history of the MW must 
be completed at early times and that no significant infall 
episode took place in the last few Gyrs, most of them based 
on the theoretical interpretation of abundances observed in 
MW disc and halo stellar populations or present-day prop- 
erties of the MW galaxy (Prantzos & Silk 1998, Boissier & 
Prantzos 1999; Hernandez et al. 2001; Alibes et al. 2001; 
Chiappini et al. 2001; Naab & Ostriker 2006; De Lucia & 
Hclmi 2008; Colavitti et al. 2008) and on kinematical argu- 
ments (e.g. Wyse 2009; Gilmore et al. 2002). 
In Fig. 9, the fact that for a few elements the predictions are 
shifted downwards with respect to the observed abundance 
ratios is mainly due to uncertainties related to the stellar 
yields used here. The [Mg/Fe] decreases by about 1.5 dex in- 
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stead of less than 0.5 dex from [Fe/H]=-2 to [Fe/H]~ 0. Con- 
cerning the elements Mg and Si, it is well known from chem- 
ical evolution studies that the yields of Woosley & Weaver 
(1995) used here tend to underestimate the [Mg/Fe]-[Fe/H] 
and [Si/Fe]-[Fe/H] relations. Francois et al. (2004) repro- 
duced the local observed abundance pattern by means of a 
model for the solar neighbourhood, increasing the yields of 
Mg by a factor of 7 and the yields of Si by a factor of 2, 
only for stars with masses m > 40Mq (see also Timmes, 
Woosley & Weaver 1995). Also the adopted C yields do not 
appear to be sufficient for reproducing the zero-point of the 
observed [C/Fe] - [Fe/H] relation and we underestimate the 
[C/Fe] on the whole metallicity range. Furthermore, The ob- 
served [C/Fe] exhibits a flat behaviour throughout the whole 
[Fe/H] range, whereas the SAM results show a remarkable 
anti-correlation. This may be certainly due to the adopted 
stellar yields, but also the IMF is likely to play some role in 
the predicted [X/Fe] - [Fe/H] relations. In Fig. 10, we show 
how different assumptions regarding the IMF and the stel- 
lar yields may affect our predictions in the [C/Fe] - [Fe/H] 
plot. With a Salpeter IMF and standard yields, our match 
to the data is very poor. The use of a Kroupa et al (1993) 
IMF causes a flattening of the predicted [C/Fe] - [Fe/H] at 
metallicities [Fe/H]> — 1. This is due to the fact that with 
the Kroupa IMF the relative fraction of intermediate mass 
stars, i.e. stars with masses 2 < M/Mq < 8 , represent- 
ing the most important C producers, is larger than with the 
Salpeter. If we assume a Kroupa et al. (1993) and we arti- 
ficially increase the yields of intermediate mass stars by a 
factor 2, all the predictions at metallicity [Fe/H]> —1 shift 
upwards and overall the predicted [C/Fe] - [Fe/H] becomes 
flatter, in better agreement with the observations. 
Romano et al. (2005) showed in detail the relative roles of 
the stellar yields and of the IMF in determining the abun- 
dance ratios, suggesting that the Salpeter IMF is not ap- 
propriate to reproduce most local chemical evolution con- 
straints. However, a detailed investigation of the stellar 
yields and of the stellar IMF requires a fine-tuning of these 
parameters and is beyond the aims of the present paper, 
which is instead focused on the chemical evolution of galax- 
ies within a cosmological framework. 

We have verified that the use of a morphological criterion 
to select MW-like galaxies has a minor impact on our re- 
sults. Among our sample of galaxies matching the chemical 
properties of MW-like galaxies, we can calculate the ratio 
between the stellar mass in the disc and the total stellar 
mass M* in the following way. The neutral gas mass of the 
MW disc is ~ 1 • 10 10 M Q (Prantzos & Silk 1998; Carroll & 
Ostlie 1996), whereas the stellar mass is ~ 5 ■W 10 M Q , with 
a ratio between the gas mass and the stellar mass ~ 0.2. We 
use this value to estimate the stellar mass in the disc of our 
selected MW-like galaxies, which is, since all the cold gas is 
in the disc, M„ jd = M g /0.2. By assuming that the stellar 
mass of the Milky Way Bulge is 1 - 2 • 1O 1O M (Ballero et 
al. 2008, Reshetnikov 2000), the observed disc-to-total stel- 
lar mass ratio is M., d /M», tot = 0.7-0.8. Of the MW-like 
galaxies selcted on the basis of their gas accretion history, 
we now select those with M* j£ j/M* jtot values ~ 0.7 — 0.8. 
In Fig. 11, we show the predicted [0/Fe]-[Fc/H] relation for 
all the MW-galaxies selected on the basis of their accretion 
history (upper panel) and on the basis of the M„^d/M* tt ot 
ratio (lower panel), and he results are nearly the same. This 



result is robust against any morphological selection criterion 
we use to select MW-like galaxies. 

In Fig. 12, we show the predicted cumulative stellar metallic- 
ity distribution (SMD hereinafter) for our MW-like galaxy 
sample, compared to observations of solar neighbourhood 
stars from various authors. The predicted SMD has been 
computed by taking into account the star formation histo- 
ries of the galaxies whose assembled mass at z ~ 2 is greater 
than 75% of their present mass. The predicted SMDs and 
the observations are normalised by requiring that the areas 
subtended by each curve have all the same values, equal to 
unity. The agreement between the observed and predicted 
SMD is remarkable. The position of the peak is much sensi- 
tive to the type la SN normalization. With our assumption 
of Aia = 0.002, the position of the peak of the SMD is satis- 
factorily reproduced. We also reproduce accurately the low- 
metallicity tail of the SMD, i.e. we predict the correct frac- 
tion of stars in the metallicity range [Fe/H]< —0.7. The high 
metallicity tail, i.e. the number o stars with [Fe/H]> 0.2, is 
in reasonable agreement with the observations. It is worth to 
note that the assumption of the IMF may have some effect 
on the predicted SMD. Fe is produced mainly by Type la 
SNe, whose rate depends on the the product of the param- 
eters k a and Ai a of eq. 6. k a depends on the IMF, hence 
a change of the IMF would require a re-tuning of the pa- 
rameter Aia in order to reproduce the correct present-day 
SN rate. However, it seems reasonable to assume that the 
quantity k a ■ Ai a is constant and fixed in order to reproduce 
the present type la SNR in MW-like galaxies. In this case, 
the IMF assumption is not relevant. 

It may be interesting to compare our results with the ones 
obtained previously by other authors. A previous relevant 
theoretical study of the chemical evolution of spiral galax- 
ies in a cosmological framework is the one of Nagashima 
& Okamoto (2006). In this paper, the authors find a good 
agreement between the predicted [O/Fc] vs [Fe/H] obtained 
for a sample of MW-like spiral galaxies and the observations 
of the abundances in local disc stars. These authors obtained 
also a good match between the predicted and the observed 
SMD. However, to model the type la SN rate and Fe produc- 
tion, Nagashima & Okamoto (2006) considered a constant 
value for the delay time of the explosion, which seems rather 
unrealistic, as indicated by a large amount of observational 
and theoretical studies of the type la SN rate in local and 
distant galaxies (Mannucci et al. 2006; Matteucci et al. 2006; 
Sullivan et al.2006; Valiante et al. 2009) which, on the other 
hand, point towards continuous delay-time distributions for 
type la SN explosions. 

Other interesting results are those obtained by Colavitti et 
al. (2008). In this paper, the authors start from cosmolog- 
ical numerical cold dark matter simulations and, from the 
merging history of the parent halos, these authors infer the 
infall history of spiral galaxies, which is implemented in a 
numerical chemical evolution code and whose evolution is 
followed a-posteriori. Colavitti et al. (2008) find an over- 
all good agreement between their predicted abundances and 
the observational constraints. However, the merging histo- 
ries are determined by means of simulations with a low time 
resolution. Furthermore, the assumption of a linear scaling 
between the baryonic accretion history and the cold-dark 
matter merging history may represent a raw approximation, 
as well as the fact that pure dark matter simulations permit 
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to explore spatial scales of the order of ~ IMpc, consider- 
ably higher than the scales of a few kpc studied in the paper 
by Colavitti et al. (2008). 

The one described in this section is the first attempt to re- 
produce the abundance ratios for such a wide set of chemical 
elements by means of a hierarchical SAM taking into account 
the lifetimes of low and intermediate mass stars and type la 
SNe. Overall, from the various results discussed in this sec- 
tion, we find that most of the chemical evolution constraints 
considered here and concerning the Milky- Way galaxy are 
satisfactorily reproduced. We find that the results discussed 
in this section are very encouraging. The agreement between 
our predictions and the observations is remarkable, in par- 
ticular concerning the [0/Fe]-[Fe/H] diagram and the SMD, 
given the very little fine tuning of the parameters performed 
so far. For other elements shown in Fig. 5, better fits to the 
data could be achievable by means of a fine-tuning of the 
adopted stellar yields, as is performed sometimes in theo- 
retical studies of the chemical evolution of the Solar Neigh- 
bourhood (e.g. Francois et al. 2004). 



3.2 Abundances in dwarf galaxies 

In this section, we present our results concerning the chem- 
ical evolution of dwarf galaxies. In particular, we focus on 
the chemical abundances in dwarf spheroidal (dSph) galax- 
ies, since for this kind of objects a large set of observational 
data is available from observations of the dSph in the Local 
Group. In this case, the selection on our galaxy catalogue 
is performed on the basis of the present stellar mass M,, 
considering all the systems with 10 6 < M*/M < 5 • 10 9 , 
broadly corresponding to the stellar mass range of local 
dwarf galaxies (Mateo 1998). In Fig. 13, we see the time 
evolution of the SFR for a few selected examples. One strik- 
ing feature of Fig. 13 is the variety of the star formation 
histories of the selected galaxies. For some systems (e.g. the 
objects labelled 18 and 20 in Fig. 13), the SFH is max- 
imum at the beginning, then smothly and constantly de- 
creasing. Some others have more than one peak (e.g. 2, 
11), mimicking the starbursty behaviour of some real dSph 
galaxies, whose complex star formation history is usually 
inferred from the observed colour-magnitude diagrams (van 
den Bergh 1994; Hernandez, Gilmore, Valls-Gabaud 2000; 
Dolphin et al. 2005). 

An important plot is the stellar metallicity vs stellar mass 
(Fig. 14), allowing us to test the star formation histories of 
our dSphs and also to tune some fundamental parameters 
regarding their chemical evolution. In this figure, the stellar 
metallicity is represented by the average stellar Fe abun- 
dance [< Fe/H >*]. This plot is of great interest since it 
is analogous to a mass-metallicity plot, with the difference 
that in general, in the mass-metallicity plots the stellar mass 
and the interstellar metallicities are represented, the latter 
usually derived by means of observations of O emission lines 
in brilliant Hji regions of star-forming galaxies (Maiolino et 
al. 2008). Local dSphs are weakly star-forming systems and 
do not allow observations of emission lines in Hu regions, 
however a plot like such of Fig. 14 provides us with simi- 
lar indications as a O/H vs M* plot, since [< Fe/H >»] is 
approximately : 



< Fe/H >*: 



J(Fe/H)(t)dM* _ J{Fe/H){t)iP{t)dt 



J dM* 



M. 



(14) 



i.e. the integral of the interstellar Fe abundance over the star 
formation history (see Thomas et al. 1999). The predicted 
[< Fe/H >»]-M» relation is compared to the one observed 
in local dwarf galaxies by Woo et al. (2008). The observa- 
tions indicate that the local dwarf galaxies clearly follow a 
mass-metallicity relation, i.e. that the smallest galaxies have 
the lowest metallicities, and that the stellar metallicity in- 
creases with stellar mass. In the left panel of Fig. 14, wc 
have assumed that dwarf galaxies can undergo mass loss 
under the same conditions as spiral galaxies, i.e. the chem- 
ical composition of the outflow is the same as the one of 
the ISM, Oi(t) = Xi(t) ■ 0{t), and for the type la SN re- 
alization probability the same value as assumed for spirals, 
Ai a = 0.002. In most of the cases, the model dwarf galax- 
ies are considerably more metal-rich than the observed ones. 
Furthermore, the predicted mass-metallicity relation is much 
flatter than the observational one. In the middle panel of 
Fig. 14, we show the predicted [< Fe/H >»]-M» relation 
computed by assuming a chemical composition of the out- 
flows as above and a very low value for the type la SN real- 
ization probability, Ai a — 0.0002, i.e. one tenth of the value 
used for MW-like galaxies. This assumption certainly low- 
ers the zero-point of the predicted mass-metallicity relation, 
however does not affect its shape, which is still flatter than 
the observed one. Both the slope and the zero-point of the 
predicted mass-metallicity relation can be further improved 
by modifying our assumptions on metal loss in dwarf galax- 
ies. In the right panel of Fig. 14, we have assumed preferen- 
tial loss of metals in the outflows, i.e. that for any chemical 
element different than H and He, the outflow rate in eq. 2 
is given by £X;(i) ■ 0(t), with £ = 10, with A Ia = 0.0002. 
As can be seen, this assumption improves considerably the 
match between the predicted and observed mass-metallicity 
relation for dwarf galaxies. One first conclusion on chem- 
ical evolution of dwarf galaxies is that, in order to repro- 
duce the slope of the observed mass-metallicity relation, one 
needs to assume that in dwarf galaxies mass loss must occur 
through metal-enhanced outflows. The importance of mass 
loss in determining the mass-metallicity relation in dwarf 
galaxies has already been ascertained (Dekel & Silk 1986). 
The fundamental role played by metal-enhanced outflows 
has already been discussed by some previous investigations 
on the mass-metallicity and luminosity-metallicity relations 
in dwarf galaxies by means of semi-analytic galaxy formation 
models (Somerville & Primack 1999). All the following re- 
sults presented in the remainder of this section are computed 
by assuming metal-enhanced outflows in dwarf galaxies. A 
more detailed discussion on the origin and importance of 
metal enhanced outflows in dwarf galaxies will be presented 
in Sect. 4. 

It may be interesting to see the effects of the star formation 
histories shown in Fig. 13 on the metallicities of the stel- 
lar populations of dwarf galaxies. In the remainder of this 
section, we will compare our predictions with observational 
data regarding local dwarf spheroidal galaxies, whose stellar 
masses are in general of the order of 10 s Mq or lower (Ma- 
teo 1998). For this reason, from this point on, in our analysis 
we will consider only the model galaxies with present stellar 
masses M* < 10 s M . In Fig. 15, we can see the predicted 
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Figure 13. Star Formation histories of a few dwarf galaxies selected according the criteria described in Sect. 3.2. 



SMDs for the same objects as those presented in Fig. 13 
. In this figure, for comparison we plot also the SMDs ob- 
served in two local dwarf spheroidals, i.e. Draco and Ursa 
Minor. The predicted SMDs have been computed with a 
number of bins comparable with the ones of the data of 
Bellazzini et al. 2002, which have 12-20 bins. For most of 
the model galaxies presented in Fig. 15, one dominant peak 
is recognizable, which corresponds to a dominant metallic- 
ity of the stellar populations of the dwarf galaxies. Some 
objects present two unresolved peaks (e.g. 6, 14), whereas 
others have two clearly distinct peaks at different metallicity 
values (11, 17). All the model galaxies with a single-peaked 
SMD or with two very close, non-resolved peaks have in 
general smooth star formation histories, or presenting small 
oscillations, or a few very short prominent peaks in their 
SFH. All the systems presenting two resolved peaks in the 
SMD have in general bursty SFHs, charactherized by quite 
broad and distinct peaks (e.g., 11, 17), or frequent short 
and high peaks overimposed to an overall continuous SFH 



(2). It is worth to note that in most of the cases, the po- 
sitions of the peaks of the predicted SMDs agree with the 
ones of the observed SMDs. It is also worth to note that 
the lowest metallicities which can be probed by our anal- 
ysis depend on the SF history. In Sect. 3.1, we have seen 
that in MW-like galaxies, we could not explore the metallic- 
ity region [Fe/H]< —2.5. MW-like galaxies have strong SF 
histories, hence a metallicity [Fe/H]~ —2.5 can be reached 
in a time shorter than the timestep used in our model. On 
the other hand, for dwarf galaxies, characterised by lower SF 
rates, the lowest metallicity probed is [Fe/H]< —3.5, reached 
sometimes on timescales higher than our timestep. Another 
reason for the low metallicities is the adopted value for Ai a , 
which in the case of dwarfs is one tenth of the value assumed 
for MW-like galaxies. 

Also in this case, owing to the large set of chemical elements 
studied here, we are able to produce predictions for the 
abundance ratios of several chemical elements. In Fig. 16, 
we show our predictions for the [O/Fe] vs [Fe/H], [Si/Fe] 
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Figure 14. Predicted distribution of the average present-day stellar metallicity in dwarf galaxies versus present stellar mass. The colour 
code, shown by the bar at the top of the Figure, represents the predicted number of galaxies with given [< Fe/H >»] with stellar mass 
Mq. The open squares with error bars are the observations in local dSphs, taken from Woo et al. (1998). In the left panel, we assume that 
the chemical composition of the ouflows Xi ou t is the same as the one of ISM and Aj a = 0.002 as assumed for MW-like galaxies. In the 
central panel, we assume chemical compositions of the outflows as above and Aj a = 0.0002. In the right panel, we assume Xi, out = £Xi, 
with £ = 10, and A Ia = 0.0002. 



vs [Fe/H] and [Mg/Fe] vs [Fe/H]. For these elements, the 
abundances are measurable in single stars of local dSphs. 
The predictions represent the abundance ratios of the stellar 
populations belonging to the selected galaxies, as described 
in Sect. 3.1. From Fig. 16, we can see how, for any ele- 
ment considered here, the characteristic decreasing trend of 
[a/Fe] vs [Fe/H] is successfully reproduced. In our models, 
this trend is due to the delay between type II SNe explosions 
and type la SN explosions, as in classical chemical evolution 
models taking into account enrichment from type la and 
type II SNe (Matteucci 2001 and references therein). 
While in the [O/Fe] vs [Fe/H] plot the observed data over- 
lap with our predictions, in the other two plots the mea- 
sured [a/Fe] are often underestimated, in particular for 
[Fe/H]> —1. The same discrepancy was found also when 



discussing the results for MW-like galaxies, and is due to the 
adopted stellar yields of Woosley and Weaver (1995), which 
lead the models to underestimate the [Mg/Fe] and [Si/Fe] 
also in the solar neighbourhood (Frangois et al. 2004). In 
these plots, a fine-tuning of the stellar yields as a function 
of the initial stellar mass would be required to reproduce the 
observations. Furthermore, we see that our models for dwarf 
galaxies predict a wider [Fe/H] range than the one spanned 
by the observational data, but this plot is not useful to un- 
derstand the relative fractions of the stellar populations of 
various metallicities. More useful in this regard is Fig. 17, 
where we show the predicted cumulative SMD of all the 
dwarf galaxies, divided in 3 stellar mass bins (see legend of 
Fig. 17), compared to the observational SMDs of the Draco 
and Ursa Minor dSph. The predicted SMDs are computed 
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Figure 15. Stellar Mctallicity distributions computed for the same set of dwarf galaxies presented in Fig. 13 (solid thick lines), compared 
with observational SMDs obtained for the Draco dSph (solid thin histograms) and for the Ursa Minor dSph (dashed thin histograms, 
sec Bcllazzini et al. 2002). 



by taking into account all the stellar populations of any se- 
lected galaxy falling in each mass bin. The cumulative stel- 
lar metallicity distributions computed for various mass bins 
are in very good agreement with the observed ones, at least 
concerning the positions of the peaks. In each mass bin, the 
predicted SMDs span very broad metallicity ranges. In par- 
ticular, in this case we predict the existence of stars at very 
low metallicities, [Fe/H]< —2.. By integrating the quantity 
dN/d[Fe/H] over different metallicity ranges, it is possible 
to estimate the relative fractions of the stellar populations 
at various metallicities. For the objects in the lowest mass 
bin, no stellar population presents [Fe/H]> —0.7. For the 
objects in the intermediate and highest mass bins, a very 
small fraction lower than 1 % have stellar populations with 
[Fe/H]> 0. 

Finally, it is interesting to note that we predict the substan- 
tial presence of very low metallicity stars with [Fe/H]< —2., 
with percentages of ~ 40%, ~ 15% and ~ 10% in the lowest, 



intermediate and highest mass bin, respectively. 
Relevant previous studies of the chemical evolution of dwarf 
spheroidals are those of Lanfranchi et al. (2006). Even if not 
carried by means of an ab-initio galaxy formation model, but 
by means of galactic chemical evolution models whose star 
formation histories are tuned from constraints of CMD ob- 
servations (see Calura, Lanfranchi & Matteucci 2008) , these 
works are useful since they allow to have some constraints 
on the parameters affecting this kind of analysis. In partic- 
ular, Lanfranchi & Matteucci (2007) present the effects of 
the wind efficiency on the predicted SMDs of dSphs, showing 
how a model where the galactic wind is suppressed leads to 
too metal-rich stellar populations. The results of Lanfranchi 
& Matteucci (2004; 2007) suggest that a strong wind effi- 
ciency is necessary to reproduce the abundance ratios and 
the stellar metallicity distributions of dwarf galaxies. 
A previous study of the chemical evolution of dwarf galaxies 
within a cosmological framework is the one of Salvadori et 
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Figure 16. Predicted distribution of the abundance ratios [X/Fe] vs [Fe/H] for dwarf spheroidal galaxies. The colour code, shown by the 
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Figure 17. Predicted cumulative stellar metallicity distribution of selected dwarf galaxies belonging to three different mass bins (solid 
thick lines) compared with observational SMDs obtained for the Draco dSph (solid histogram) and for the Ursa Minor dSph (dashed 
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al. (2007) . These authors model chemical evolution of dwarf 
galaxies by taking into account the stellar lifetimes of low 
and intermediate mass stars and of massive stars, however 
they do not consider type la SN explosions. Salvadori et 
al. (2007) explain the observed decline of [a/Fe] vs [Fe/H] 
in local dSphs by means of differential galactic winds. This 
result is apparently in contrast with a large variety of chemi- 
cal evolution studies, which underline the importance of the 
time delay of type la SN explosions in explaining the de- 
crease of the [a/Fe] vs [Fe/H] in local dSphs (Recchi et al. 
2001; Ikuta & Arimoto 2002; Lanfranchife Matteucci 2003; 
Lanfranchi et al. 2006; ; Recchi et al. 2006; Marcolini et al. 
2008). 

Very recently, Sawala et al. (2009) studied the chemical evo- 
lution of local dwarfs in a cosmological framework by means 
of high-resolution hydrodynamical numerical simulation, in- 
cluding also type la and type II SN enrichment. Their results 
indicate that the chemical properties of dwarf galaxies are 
driven by several aspects, such as the SN feedback, the UV 
background and their gravitational potential. Their results 
stress the importance of efficient metal-enhanced outflows, 
fundamental for reproducing the stellar mass-metallicity re- 
lation. 

A numerical study of the stellar metallicity distributions of 
dwarf spheroidals has been performed by Ripamonti et al. 
(2007). Their results point towards the substantial presence 
of very low metallicity [Fe/H]<-3 stars in dwarf spheroidals, 
in agreement with the results found in this paper. If no ob- 
servational bias against these stars is present, possible so- 
lutions to this problem invoke a primordial IMF truncated 
for stellar masses below 1 Mq or a pre-enrichment of the 
gas. However, very recently, the presence of extremely low 
metallicity stars has been detected in ultrafaint dwarf galax- 
ies (Frebel et al. 2009). 

3.3 Abundances in elliptical galaxies 

We conclude this section presenting our chemical evolution 
results obtained for elliptical galaxies. Elliptical galaxies are 
selected on the basis of their present-day (B-V) colour. In 
specific, we consider elliptical galaxies all the systems with 
present (B-V)> 0.85, following Roberts & Haynes (1994). 
We perform a further cut on the basis of the present- 
day velocity dispersion a, considering all the systems with 
30 < log(a/km/s). This range for the velocity dispersion 
is similar to the one spanned by the observational data of 
Thomas et al. (2005), who determined observationally the 
integrated stellar chemical abundances in a sample of local 
field ellipticals, and whose set of data will be used here for 
comparison with our predictions. In Fig. 18, we show the 
star formation histories of a few ellipticals selected accord- 
ing to the criteria described above. The numbers present in 
each panel of Fig. 18 indicate the present-day stellar masses 
of the selected galaxies. The selected SFRs are characterised 
by numerous peaks occurring in most of the cases at cosmic 
times < 6 Gyr, corresponding to redshift z > 1. In general, 
galaxies with larger stellar masses exhibit higher SFR val- 
ues. The SFHs of the objects shown in Fig. 18, with little or 
no substantial SF in the last few Gyrs, are qualitatively in 
agreement with the SFHs of early type galaxies drawn from 
other SAMs (De Lucia et al. 2006; Pipino et al. 2009). 
Generally, in local ellipticals the stellar abundances are mea- 



sured by means of absorption- line indices, such as Mg6 and 
< Fe >= 0.5(Fe52720 + Fe5335) (see Thomas et al. 2005 
and references therein). This means that the observational 
abundances are luminosity-averaged values, which in prin- 
ciple represent underestimates to the true, mass-averaged 
abundances, which are the ones computed by means of our 
models, because of the fact that at constant age, metal-poor 
stars are brighter (Greggio 1997). However, as discussed in 
several papers (Matteucci, Ponzone & Gibson 1998; Thomas 
et al. 1999; Recchi et al. 2009), detailed chemical evolution 
calculations show that the difference between mass-averaged 
and luminosity-averaged abundances is very small, typically 
of the order of a few 0.01 dex. This result is not dependent 
on fundamental chemical evolution parameters, such as the 
stellar IMF (Recchi et al. 2009). For this reason, we neglect 
this difference and we compare our predicted mass-averaged 
stellar abundances to the ones derived observationally by 
means of absorption-line indices. 

In the upper panel of Fig. 19, we show the predicted 
average stellar a/Fe ratio versus velocity dispersion, com- 
puted assuming a Salpeter IMF and compared with a set 
of local observations by Thomas et al. (2005). In this case, 
the predicted [< a/Fe >*] is computed considering O as 
representative of a— elements. To compute our predictions, 
in Fig. 19 we use a standard Salpeter IMF, and a type 
la SN realization probability Ai a = 0.002. As we will see 
later, with our standard assumptions, this value allows us 
to correctly reproduce the zero point of the stellar [a/Fe] - 
a relation. 

While the observations indicate a positive correlation be- 
tween [a/Fe] and a, with galaxies with larger a presenting 
larger a— enhancement, our predictions indicate a flat [a/Fe] 
-a relation. In general, the observed [a/Fe] -a is interpreted 
in terms of smaller formation timescales in larger ellipti- 
cals, in which the chemical enrichment is dominated by the 
contribution of type II SNe, hence most of the stars form 
with a chemical composition rich of a elements (Matteucci 
1994; Thomas et al. 2005). On the other hand, smaller galax- 
ies present lower stellar [a/Fe] since they form on longer 
timescales, comparable to the exploding times of type la 
SNe, whose contribution to the Fe enrichment of the ISM 
is relevant, lowering the average stellar [a/Fe]. Matteucci 
(1994) has shown that a monolithical galaxy formation sce- 
nario for ellipticals is successful in reproducing the observed 
stellar abundance ratios. According to this scenario, in el- 
lipticals star formation stops owing to the onset of galactic 
winds, which develop on smaller timescales in larger ellip- 
ticals, according to an inverse-wind scheme. Calura et al. 
(2008) have confirmed that gas outflows should occur on 
smaller timescales in large galaxies, according to a "down- 
sizing" pattern for galaxy formation. 

The impossibility to reproduce the observed [a/Fe] -a re- 
lation is a well-known problem, common to many hierar- 
chical SAM. The reason is that, in hierarchical models, the 
star formation of large elliptical galaxies is extended until 
relatively recent times, hence without avoiding that large 
amounts of Fe are restored by type la SNe into the ISM 
while star formation is still active. In general, this leads 
large galaxies to present lower average stellar [a/Fe], i.e. 
to an anti-correlation between [a/Fe] and a (Thomas 1999; 
Nagashima et al. 2005; Pipino et al. 2009), contrary to what 
the observations indicate. On the other hand, by assuming 
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Figure 18. Star Formation histories of a few elliptical selected according the criteria described in Sect. 3.3. 



a simple Salpeter IMF constant in time, our results indicate 
rather a flat relation between these two quantities. This is 
still at variance with the observations, however, this should 
be certainly regarded as a an improvement with respect to 
the previous attempts discussed above. 

Our basic assuptions, i.e. a Salpeter IMF, the inclusion of 
interaction-triggered starbursts and AGN feedback do not 
allow us to produce a positive correlation between the stel- 
lar [a/Fe] and a. In a recent paper, Arrigoni et al. (2009) 
studied the chemical evolution of local ellipticals by means 
of a semi-analytic model for galaxy formation which includes 
AGN feedback. Their main result concerns the capability to 
reproduce the observed correlation between stellar [a/Fe] 
and stellar mass, after a fine-tuning of the a type la SN con- 
stant and by assuming an IMF slope x = 1.1, where the 
canonical Salpeter slope is x — 1.35. As already stressed, in 
chemical evolution studies the type la SN constant Ai a is an 
uncertain parameter, however the IMF is in general better 
constrained by local observations. In the light of the results, 
now we aim at studying the [a/Fe] - a relation in the same 
conditions as the ones of Arrigoni et al. (2009), i.e. assum- 
ing a flatter IMF. In the lower panel of Fig. 19, we show 
the stellar [a/FcJ-a relation predicted for massive galaxies, 
obtained by considering a constant intial mass function with 
index x = 1. This figure shows clearly that the assumption 
of a constant IMF flatter than the Salpeter has the effect of 
lifting the zero-point of the predicted [a/Fe] -a relation, but 
has no effect on its slope. This is also in agreement with the 
discussion about the role of the IMF in the chemical evolu- 
tion of early type galaxies of Pipino et al. (2009). 
Recent observational and theorethical results seem to in- 
dicate a dependency between the star-formation rate and 



the slope and the upper mass limit of the IMF. In normal 
star formation conditions, i.e. in local regions with SFRs 
< lOOMo/j/r, a standard IMF slope x ~ 1.3 for stellar 
masses > IMq is commonly accepted (Kroupa 2002, Rec- 
chi et al. 2009). On the other hand, in local starbursts, in 
the cores of stellar clusters and in ultra-compact galaxies, 
with SFRs > 100Mo/yr, various results indicate a flatter 
IMF, with a slope x < 1 (Elmegreen 2009; Dabringhausen 
et al. 2009). Further evidences in favour of a flat IMF in 
environments undergoing strong star formation come from 
chemical evolution models, which, in order to be able to 
explain the oversolar [a/Fe] values observed in the hot in- 
tracluster gas, must necessairily invoke an IMF skewed to- 
wards massive stars (Gibson & Matteucci 1997, Portinari et 
al. 2004). Motivated by these results, now we aim at inves- 
tigating the impact of a SF-dependent IMF on the stellar 
[a/Fe] vs a relation predicted by local ellipticals. Our aim 
is also to test how our results depend on other ingredients 
of the SAM model used in this paper, such as interaction- 
triggered starburst in massive galaxies and AGN feedback, 
and the importance of including or neglecting these effects 
in conjuction with the main chemical evolution parameter, 
i.e. the IMF. 

In the left panel of Fig. 20, we show the predicted [a/Fe]- 
<j relation obtained by considering a SF-dependent IMF, 
characterized by a standard Salpeter slope (x = 1.35) in 
objects with SFRs < 100M Q /yr, and x = 1 in objects with 
SFRs > 100Mo/yr. This assumption in consistent with the 
study of Dabringhausen et al. (2009), indicating IMF slopes 
x < 0.7 in starburst and ultra-compact galaxies, with SFRs 
of 10 — lOOM0/yr. Here we consider an IMF slope x = 1, 
close to the upper limit of Dabringhausen et al. (2009), since 
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Figure 19. Predicted stellar average a/Fe ratio vs velocity dispersion compared to local observations by Thomas et al. (2005) (open 
diamonds) . The colour code, shown by the bar at the top of the Figure, represents the predicted number of galaxies with given [< a/ Fe > * ] 
with velocity dispersion cr, normalised to the total number of galaxies with that a. In the upper panel, we show our results computed 
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Figure 20. Predicted stellar average a/Fe ratio vs velocity dispersion (contour), compared to local observations by Thomas et al. (2005) 
(open diamonds). The colour code represents the predictions, as described in Fig. 19. In the left panel, our results have been computed 
by means of a SF-dependent IMF and our standard assumptions, which include starburst-triggered interactions at high redshifts and 
AGN feedback. In the right panel, we have assumed a a SF-dependent IMF and AGN feedback but we have supressed starburst-triggered 
interactions. 



this does not require substantial modifications of the pa- 
rameters of our SAM to reproduce the main local observa- 
tional constraints, such as the luminosity function or the 
Tully-Fisher relation. In this case, we see that the use of 
this IMF allows us to considerably improve our results, re- 
producing the observed correlation between [a/Fe]-cr and 
bringing our results in good agreement with the observa- 
tions on a large range of velocity dispersions. In the right 
panel of Fig. 20, we show our results obtained with a model 
which does not include encounter-triggered starbursts, but 
it does include a SF-varying IMF as described above. In this 
case, the predicted [a/Fe]-er relation is too shallow to repro- 
duce the observations. This result outlines the importance 
of interaction-driven starbursts in determining the [a/Fe]-a 
relation of local elipticals. At this stage, it may be interest- 
ing to test also the role of the AGN feedback in the chemical 



evolution of ellipticals. In Fig. 21, we compare our results ob- 
tained with our standard model and with a SF-dependent 
IMF with the ones achieved with a model which does in- 
clude interaction-triggered starbursts, but which does not 
take into account the effect of AGN feedback in massive 
galaxies. The model with no AGN feedback has a shallower 
slope than the standard model and produces a poor match 
to the observations. This indicates that also the AGN feed- 
back must play some important role in shaping the [a/Fe]-a 
relation of local ellipticals. 

In summary, in order to reproduce the correlation between 
the stellar [a/Fe] and a of local ellipticals with our SAM, 
three basic ingredients seem to be required: (i) encounter- 
triggered starbursts at high redshift; (ii) AGN feedback; (iii) 
a SF-dependent IMF, flatter than the Salpeter in systems 
with strong SF. It is worth to note how a small decrease of 
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Figure 21. Predicted stellar average a/Fe ratio vs velocity dispersion (contour), compared to local observations by Thomas et al. (2005) 
(open diamonds). The colour code represents the predictions, as described in Fig. 19. In the left panel, our results have been computed 
by means of a SF-dependent IMF and our standard assumptions, which include starburst-triggered interactions at high redshifts and 
AGN feedback. In the right panel, we have assumed a a SF-dependent IMF and starburst-triggered interactions but we have supressed 
the effect of AGN feedback. 



the IMF slope allows us to considerably improve our results. 
A further decrease of the IMF index x may possibly improve 
our results, but, in order to reproduce the local constraints, 
it would probably require a re-calibration of the main pa- 
rameters of the SAM to reproduce the local constraints. This 
will be the subject of future work. 

The predicted mass-stellar metallicity relation is shown in 
Fig. 22, compared to the same relation as observed in lo- 
cal ellipticals (Thomas et al. 2005). In the upper panel of 
Fig. 22, we present our results obtained with a Salpeter 
IMF. In this case, we predict a basically flat [< Z/H >*] 
vs a relation, whereas the data indicate a positive correla- 
tion between these two quantities. In the lower panel of Fig. 
22, the [< Z/H >»] -a relation is computed by assuming a 
SF-dependent IMF as described above. Again, this assump- 
tion improves considerably our fit to the observational data. 



This plot povides a further, independent indication that the 
assumption of a SF-dependent IMF is a possible solution to 
the long-standing problem of explaining the chemical abun- 
dances in ellipticals with hierarchical SAMs. 
When comparing our predictions with observations, it is 
worth to stress that the observational [< Z/H >»] vs a 
is likely affected by a template bias (Proctor et al. 2004; 
Thomas et al. 2005; Nagashima et al. 2005), which reflects 
the fact that the stellar population models used to determine 
the metallicities are based on stellar spectra in which some 
stars present a sub-solar [a/Fe]. This effect could cause a 
steep [< Z/H >*] vs a relation, but, as shown by Thomas et 
al. (2005), it leaves basically unchanged the [< a/Fe >»] vs 
a relation. Furthermore, as outlined by Pipino et al. (2009), 
the fact that elliptical galaxies exhibit strong metallicity gra- 
dients makes a detailed comparison between our predicted 
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[< Z/H >»] vs a relation with observations quite difficult, 
and is likely to mimic a steeper [< Z/H >»] vs a relation, 
leaving essentially unchanged the [< a/Fe >»] vs a. 
Finally, it is important to stress that in this paper, to com- 
pute the integrated stellar abundances in early-type galax- 
ies, we have considered the integral SF histories of our se- 
lected galaxies, represented by the sum at each timestep 
of the SF histories of all the progenitors. As stressed by 
Pipino et al. (2009), the most massive progenitors have in 
general shorter star formation timescales and should dom- 
inate the integrated stellar abundances. For some systems, 
this may cause an underestimation of the stellar [a/Fe] ra- 
tios by 0.1-0.2 dex. This effect is not taken into account at 
the present time, as considering in detail the star formation 
histories of the single progenitors for all our galaxies would 
require considerable computational times. A subsequent pa- 
per completely dedicated to this aspect is currently under 
preparation. 



4 DISCUSSION AND CONCLUSIONS 

The major aim of the present paper is to model chem- 
ical evolution by means of a hierarchical SAM by taking 
into account properly all the various stellar sources. We also 
aim at providing chemical evolution predictions computed 
by means of star formation and infall histories derived from 
first principles, by means of an ab-initio approach. Beside 
high-mass stars, dying as core-collapse type II SNe, we con- 
sider also the contributions from type la SNe and low and 
intermediate mass stars. A continuous delay-time distribu- 
tion is assumed for type la SNe, relevant producers of some 
important elements, such as Fe and Si. We also consider fi- 
nite lifetimes for low and intermediate mass stars, which are 
important producers of C and N. Our methods allow us to 
provide, for the first time, chemical evolution predictions for 
a large set of chemical elements, produced by stars of vari- 
ous types on different timescales. This is particularly useful 
to study the abundance ratios between elements synthesised 
by different stellar sources on various timescales, which are 
helpful tools to measure the timescales of star formation 
in galaxies of various masses and morphologies. We stud- 
ied chemical evolution in various galactic types, comparing 
our predictions with available observations of chemical abun- 
dances in the Milky Way, in local dwarf galaxies and in local 
ellipticals. 

First, we studied the abundance ratios vs metallicity for 
Milky Way-like galaxies. We compared our predictions to 
observations available for local stars of the Solar Neighbour- 
hood for various elements. Overall, concerning the most im- 
portant constraints, we found a good agreement between our 
predictions and the observations. In the [O/Fe] vs [Fe/H] 
plot, our results for MW-like galaxies indicate the presence 
of a horizontal population of stars with metallicity —2 < 
[Fe/H] < -0.2 and constant [0/Fc]~ -0.2- -0.3. This pop- 
ulation is the result of substantial late infall episodes of 
pristine gas, very frequent in the selected MW-like galax- 
ies, which have the effect of decreasing the metallicity of the 
most recently formed stellar populations, leaving unchanged 
their [O/Fe] ratio. When we consider only the galaxies with 
at least 75 % of the present-day mass already assembled 
at z = 2, this horizontal population of stars with constant 



[0/Fe]~ —0.3 disappears. These galaxies do not experience 
major accretion episodes at epochs after ~ 6 Gyr. If we con- 
sider only these subset of galaxies, the observational [O/Fc]- 
[Fe/H] is accounted for. Our results concerning the [O/Fe] 
vs [Fe/H] plots exclude that our Galaxy experienced major 
gas accretion episodes within the last few Gyrs, either due 
to merging with gas-rich, low-metallicity dwarf galaxies, or 
due to infall of pristine gas. This is is in agreement with a 
large set of previous results based on the interpretation of 
abundances observed in MW disc and halo stellar popula- 
tions and on kinemathical studies of the various components 
of the MW (e.g. Gilmore et al. 2002; Colavitti et al. 2008 
and references therein). 

The typical decrease of the observed [a/Fe] vs [Fe/H] plot 
is reproduced, as well as the observed trends for the other 
abundance ratios vs [Fe/H]. However, for other elements, 
the predictions are shifted downwards with respect to the 
observed abundance ratios. This is due to the adopted stel- 
lar yields, which in most cases lead to underestimating the 
observed abundance pattern. For elements other than O, 
the uncertainty in the stellar yield is a well-known prob- 
lem of chemical evolution models of the solar neighbourhood 
(Frangois et al. 2004). Also the IMF is likely to play some 
effects on the [X/Fe] vs [Fe/H] plots. For instance, with a 
Salpeter IMF, we predict an anti-corrlation between [C/Fe] 
and [Fe/H], whereas the observational data indicate a sub- 
stantially flat behaviour. We have shown how, by assuming 
the IMF of Kroupa et al. (1993) and by slightly modifying 
the stellar yields of C for both intermediate mass- and mas- 
sive stars, it is possible to improve our fit to the [C/Fe] - 
[Fe] diagram. From the study of the local stellar metallicity 
distribution, we showed that we can correctly reproduce the 
position of the peak and the low-metallicity tail of the ob- 
served SMD. Also the predicted fraction of stars with metal- 
licities < [Fe/H] < +0.2 is in reasonable agreement with 
the observational data. 

By assuming that the chemical composition of the outflows 
is the same as the one of the ISM and the same value for 
type la SN realization probability Ai a as assumed for MW- 
like galaxies, we severely overestimate the stellar metallic- 
ities of dwarf galaxies. Differential outflows, where metals 
are lost more efficiently than H and He, have a major effect 
on the slope of the predicted < [Fe/H] >„ — M* relation, 
whereas the quantity Ai a affects primarily the zero point 
of this relation. The assumptions of strongly enhanced out- 
flows and of a lower value for Ai a substantially decrease the 
the average stellar metallicity of dwarf galaxies, yielding a 
better match between predictions and observations in the 
< [Fe/H] >* — M» diagram. Other theoretical results out- 
line the role of metal-enhanced outflows in dwarf galaxies 
(see Recchi et al. 2008). The physical justification for differ- 
ential mass loss from dwarf galaxies has been studied in the 
past by means of numerical simulations. MacLow & Ferrara 
(1999) modelled the effects of SN explosions in dwarf galax- 
ies, finding that the metal ejection is much more efficient in 
lower mass galaxies. By means of 2D hydrodynamical sim- 
ulations, D'Ercole & Brighenti (1999) studied the feedback 
of a starburst on the ISM of typical gas-rich dwarf galaxies. 
They showed that metals are expelled more easily than the 
global ISM since the metal-rich material shed by the massive 
stars belongs to a very hot phase of the ISM, which can be 
easily accelerated to velocities higher than the escape speed 
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Figure 22. Predicted stellar average Z/H vs velocity dispersion <r, compared to local observations by Thomas et al. (2005) (open 
diamonds). In the upper and lower panel, we show the results obtained by assuming a Salpeter IMF and a SF-dependent IMF (i.e. with 
x = 1.35 for systems with SFR< IOOMq /yr and with x = 1 for systems with SFR< IOOMq /yr), respectively. The colour code, shown 
by the bar at the top of the Figure, represents the predicted number of galaxies with given [< Z/H >*] with velocity dispersion a, 
normalised to the total number of galaxies at the considered a. 



and leave the galaxy. On the other hand, the interstellar gas, 
heated up by the SN explosions is only temporarily affected 
by the starburst, and the galaxy is able to recover a cold ISM 
after a time of the order of ~ 100 Myr from the starburst. 
Also results from chemical evolution models suggest that 
heavy elements should be lost more effectively than H to 
reproduce the global properties of dwarf galaxies (Pilyugin 
1993, Marconi et al. 1994). The existence of metal-enhanced 
outflows from star-forming galaxies has also been recently 
confirmed by observations of local starbursts (Martin et al. 
2002; Ott et al. 2005). 

Our results concerning the abundance ratios of dwarf galax- 
ies are compared to the abundances observed in local dSphs. 
Also in this case, the characteristic decreasing trend of 
[a/Fe] vs [Fe/H] is reproduced, and explained as due to the 
delay between type II SNe explosions and type la SN explo- 
sions. 

The study of the individual and cumulative stellar metallic- 
ity distributions of dwarf galaxies of various masses provides 
a satisfactory agreement between predictions and observa- 
tions. Also this result is a direct consequence of having as- 
sumed metal enhanced outflows in dwarf galaxies and a low 
value for Ai a . The only discrepancy concerns the fraction 
of stars with metallicity [Fe/H]< —2. Our predictions in- 
dicate that they must be present in substantial fractions, 
in particular in galaxies with the lowest masses. The ob- 
servational SMDs used in this paper do not indicate the 
presence of these low metallicity stars in dwarf spheroidals. 
However, very recently a few extremely low metallicity stars 



have been detected in the Sextans and Draco dSphs and in 
ultrafaint dwarf galaxies (Aoki et al. 2009; Cohen & Huang 
2009; Frebel et al. 2009). 

Finally, we studied the chemical evolution of local ellipti- 
cals. Our results were compared with the available observa- 
tions, which represent the average stellar abundances in local 
early-type galaxies. Concerning the [a/Fe] vs a plot, the ob- 
servations indicate a positive correlation between [a/Fe] and 
a, usually interpreted in terms of more efficient star forma- 
tion efficiency in larger galaxies (Matteucci 1994; Thomas et 
al. 2005). Previous attempts to reproduce the observations 
in early type galaxies by means of semi-analytical ACDM 
galaxy formation models have led to several insuccesses. 
Theoretical results by various authors point towards an anti- 
correlation between a/Fe and a (Thomas 1999, Nagashima 
et al. 2005). 

We have performed a detailed study of the role of various 
ingredients in determining the slope and the zero point of 
the a/Fe - a relation. Our predictions computed with a con- 
stant Salpeter IMF indicate a flat [a/Fe] -a. We have found 
that, by assuming a SF-dependent IMF, Salpeter-like (i.e. 
with slope x = 1.35) in objects with SFR < IOOMq /yr and 
slightly flatter (with x = 1) in systems with strong star 
formation, the observed correlation can be accounted for, at 
least on the velocity dispersion range 1.8 < log (a / 'kms / s) < 
2.2. A major role is played by starburst-triggered interac- 
tions and AGN feedback. In fact, we have seen that when 
we relax any of these two ingredients in our models, we pre- 
dict a flatter [a/Fe] -a relation than the one produced con- 
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sidering both effects and the SF-dependent IMF. Our main 
result is then that, in order to reproduce this correlation, 
all these three ingredients must be taken into account. The 
discrepancy between our predicted [a/Fe] -a relation and 
the observations for log(a /kms/s) > 2.3 may be due to the 
fact that, at the present stage, the integrated stellar [a/Fe] 
have been computed by considering the integral SF histo- 
ries of early type galaxies, and not the contributions of the 
single progenitors. Since in general the most massive pro- 
genitors should exhibit the highest stellar [a/Fe] and have a 
higher weight in the calculation of the total [a/Fe], a more 
detailed computation of the stellar [a /Fe] values taking into 
account each single progenitor could further steepen the pre- 
dicted [a/Fe]-cr relation and alleviate the discreancy for the 
largest mass systems. This task requires long computational 
times and will be the subject of a forthcoming paper. With 
a standard Salpeter IMF, we cannot reproduce the slope of 
the mass-stellar metallicity relation observed in ellipticals, 
which should however be considered a less robust observa- 
tional constraint, mainly owing to a possible template bias 
affecting the stellar population models used to determine 
the metallicities. However, we have shown that, also in this 
case, the assumption of a SF-dependent IMF alleviates the 
discrepancy between our predictions and the observational 
data. 

In the future, our semi-analytical model for galaxy forma- 
tion will be used to perform a more thorough study of the 
abundance ratios in ellipticals and of the mass-metallicity 
relation and its evolution (Calura et al. 2009). 
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